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We show that a model of interacting electrons in one dimension is able to explain the order of 
magnitude as well as the temperature dependence of the critical supercurrents recently measured 
in nanotube samples placed between superconducting contacts. We use bosonization methods to 
compute superconducting correlations in the presence of the long-range Coulomb interaction, ending 
up with a picture in which the critical current does not follow the temperature dependence of the 
gap in the contacts, in contrast to the prediction of the conventional proximity effect. Our results 
also reveal the presence of a short-range attractive interaction in the nanotubes, which accounts for 
a significant enhancement of the critical supercurrents. 
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Since the discovery of carbon nanotubes, these have of- 
fered a great potential for novel electronic properties and 
technological applications. It has been checked experi- 
mentally Q the theoretical prediction that there should 
be semiconducting as well as metallic nanotubes 0. It 
has been also remarkable the experimental observation of 
unconventional transport properties j| , that seem to be 
compatible with the expected Luttinger liquid behavior 
of one-dimensional electron systems Q]. Different ap- 
proaches have predicted the appearance of phases with 
broken symmetry in the carbon nanotubes at very low 
energies []|-|| . Anyhow, the estimates are in general that 
there should be enough margin to observe the character- 
istic scaling behavior of the Luttinger liquid over a wide 
range of temperatures. 

A different class of experiments has been aimed to test 
the superconducting properties of the carbon nanotubes 
P,^0[. One of the most striking results has been the ob- 
servation of supercurrents along carbon nanotubes placed 
between superconducting contacts || . In a sample made 
of a single-walled nanotube, for instance, critical super- 
currents have been measured that are about 40 times 
higher than expected from the value of the gap in the 
contacts || . They also show a very flat dependence with 
temperature, until the critical value of the superconduct- 
ing contacts is approached. In that respect, there is a 
marked difference from the behavior of another sample 
made of a rope of nanotubes, where the critical super- 
current seems to follow the BCS gap in a certain range 
of low temperature Q . 

A model of the electron interaction in the carbon nan- 
otubes should give a quantitative account of all these dif- 
ferent observations of superconducting correlations. We 
show in this letter that the mentioned features of the crit- 
ical supercurrent can be understood in the framework of 
a one-dimensional theory of interacting electrons. We 
will see that the experimental data are consistent with a 
definite form of the one-dimensional interaction, as it is 



actually nontrivial to reproduce both the shape and the 
order of magnitude of the supercurrents in the single- 
walled nanotube and in the rope of nanotubes. In par- 
ticular, the experimental values of the supercurrent point 
at a sensible renormalization of the strength of the long- 
range Coulomb interaction, specially in the sample made 
of a rope of nanotubes, in agreement with earlier theo- 
retical predictions Jll],[l2j • 

A metallic single-walled nanotube has several one- 
dimensional subbands, with two pairs of linear branches 
crossing at Fermi points hp and — hp. We deal in this 
letter with an effective description of the nanotubes for 
energies below the scale E c at which all the gapped sub- 
bands decouple in the computation of low-energy proper- 
ties, so that the relevant modes left belong to the linear 
branches close to the Fermi level. We can estimate this 
energy E c as a few tenths of eV, for a typical single- walled 
nanotube with about 10 subbands. 

The low-energy excitations can be encoded into four 
boson fields, each boson corresponding to a linear branch 
in the same fashion as in the Luttinger model Q. The 
hamiltonian of the effective theory can be written in 
terms of the respective density operators pi a <n labelled 
by the Fermi point a = 1, 2 and by the chirality i — L, i?, 
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In the above expression, k c is related to E c through the 
Fermi velocity vp, k c = E c /vp. 

Our assumption regarding the interaction will be the 
presence of the long-range Coulomb interaction V(k) « 
e 2 /(47r 2 ) log |fc c /fc|, which remains unscreened in one spa- 
tial dimension |l3|], plus an additional short-range effec- 
tive attraction coming from the coupling to the elas- 
tic modes of the nanotube. In this framework, we are 
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neglecting backscattering and Umklapp processes that 
mix different chiralities and Fermi points, relying on the 
fact that those interactions have smaller relative strength 
(~ 0.1 /n, in terms of the number n of subbands 14 1) 
and they remain small down to extremely low energies 

' 

The correlators in the model governed by (m can be 
computed by changing variables to the total charge den- 
sity operators 



i = L,R 



(2) 



where N stands in general for the number of channels 
{(a, a)}, so that N = 4 in the case of a single- walled 
nanotube. 

A typical propagator of Cooper pairs, for instance, be- 
comes 



G(x,t) = (* Lir (x,t)* ml (x,t)*ti r (0,0)*t2iM) 



C(x,t)F(x,t) 



(3) 



where F is the part that does not depend on the inter- 
action and C corresponds to the propagation of the total 
charge. At zero temperature, for instance, we have 



C(x,t) 
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where fi(k) = 1/^1 + 2NV(k)/v F and v F 
The other factor has the simple dependence 

N- 

F{x,t) = 1/ \k*(x- v F t)(x + v F t)\ N 



(1 — cos(fcx) cos(v F kt)) 
(4) 

v F /n{k). 
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The critical supercurrent X can be estimated under the 
assumption that (i) the normal-superconductor junctions 
are perfectly transmitting |]l5| , or (ii) the single-particle 
scattering is relevant at the interfaces |[6]]. The latter 
is more realistic for the experiments that we are con- 
sidering. The distance L between the superconducting 
contacts is large enough that X can be expressed as a 
function of L and the temperature T as 

,1/T 

X L {T) = ev F k c dTG{L,-ir) (6) 
Jo 

In the above equation, G stands for the appropriate ex- 
pression at finite temperature. The analytic continuation 
to imaginary time, however, cannot be taken directly in 
expressions like (Q), and for computational purposes it 
is more convenient to introduce the temperature depen- 
dence through the Matsubara formalism 



C(x, — ir) — exp 
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1 — cos(fca;) cos(27rrnTr) 
(27rmT/w F ) 2 + k 2 
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We can use Eq. (||) to test whether the behav- 
ior of the critical currents measured in Ref. || can be 
reproduced within the present framework. The com- 
parison should be fairly direct for the sample that is 
made at one end of a single nanotube (called ST\ 
in Ref. According to the above discussion, we 

consider a momentum-dependent parameter fi(k) — 
+ N (e 2 /(2ir 2 v F ) log \k c /k\ - g/(-KV F )) 7 taking in 
this case a number of channels N — 4 in the above equa- 
tions. 

We have checked first that the critical current X s» 
0.1 /i A of the STi sample atTw OK, found anomalously 
high in the BCS framework, can be explained with the 
present model. The results represented in Fig. [j] show 
the magnitude of Xl (0) for different values of the inter- 
action, the distance being measured in units of k~ l . The 
actual values of the supercurrent are obtained by multi- 
plying the magnitudes in Fig. [I] by the prefactor in Eq. 
(^). The Fermi velocity can be obtained from the hop- 
ping amplitude t i=s 2.1eV and the nearest-neighbor dis- 
tance a s» 1.4 A, by using the expression v F = 3ta/2. A 
reasonable estimate of the cutoff k c for the single- walled 
nanotube is k c 0.5nm , which gives ev F k c 30/iA. 




FIG. 1. Plots of the critical current (in units of ev F k c ) ver- 
sus distance, at T — 0, for different strengths of the Coulomb 
interaction. From top to bottom, the solid curves correspond 
to 2e 2 /(ir 2 v F ) = 1.0,2.0,4.0,8.0. The dotted (colored) lines 
correspond in each case to the correction by effect of the ad- 
ditional short-range interaction, with 4g/(irv F ) = 0.75. 

We observe that the coupling corresponding to the 
bare parameters of a graphite layer, 2e 2 /(tt 2 v f ) w 8.0, 
does not lead to sensible results. The total length of the 
sample STi is ~ 300nm, and we should expect to get 
the correct order of magnitude of the critical current at 
L ~ 50/fc c . The most appropriate value for the Coulomb 
interaction seems to be 2e 2 / '{tt 2 v f ) w 1.0. The order 
of magnitude X w 0.1/LtA is then reached, most precisely 
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if one takes into account a coupling for the short-range 
attractive interaction gj (ttvf) ~ 0.2. The sensible reduc- 
tion in the value of 2e 2 / (tt 2 vf) can be understood by the 
presence of nearby charges and the renormalization of the 
interaction in the narrow rope into which the nanotubc 
merges, as we discuss afterwards. 

Moreover, the dependence of the critical current on 
T for the mentioned interaction strength reproduces the 
shape that has been observed in the measurements of the 
sample ST\. In the theoretical model, the temperature 
T is given in units of the unique energy scale E c . This 
means that the critical temperature T c pa 0.4K of the con- 
tacts for the sample ST\ corresponds to the dimensionless 
value T c /E c pa 2 x 10~ 4 . We have represented in Fig. || 
the results for the critical current Tl(T) at L — 50/fc c , 
with and without the effect of the short-range attractive 
interaction. It is remarkable the smooth behavior of the 
critical current in the low-temperature regime below T c . 
The slight increase observed near zero temperature is re- 
lated to the renormalization of the transmission at the 
interfaces fl7}| . Near T c , the suppression of superconduc- 
tivity in the contacts should be incorporated to produce 
a sharp decrease, leading then to the full agreement with 
the experimental results of Ref. 0] . 
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FIG. 2. Plots of the critical current (in units of evFk c ) 
versus T/E c , for 2e 2 /(tt' 2 vf) — 1.0 and a coupling of the 
short-range attractive interaction 4g/(iTVF) = (lower curve) 
and 0.75 (upper curve). 

Moving now to the sample made of a rope of nanotubes 
(called RO3 in Ref. H), we have to take into account two 
main differences with respect to the preceding discussion. 
First, the energy scale E c up to which the rope can be 
seen as a purely one-dimensional system is smaller, com- 
pared to the cutoff introduced for the sample ST\ . Given 
that the diameter of -RO3 can be approximately 15 times 
larger than that of the single- walled nanotube of STi , we 
may assume that the energy cutoffs in the two samples 



differ in the inverse proportion by a factor of 15. This 
means that a given temperature of the sample RO3, when 
measured in units of the corresponding E c , looks compar- 
atively higher than the same temperature in the sample 
ST\. Thus, the critical temperature T c pa 1.1K of the 
contacts used for the sample .RO3 gives a dimensionless 
value T c /E c pa 9 x 1CT 3 , which is more than one order of 
magnitude higher than the ratio for the sample ST%. 

The second important difference between the rope of 
nanotubes RO3 and the sample ST± is the interaction 
among the large number of nanotubes (pa 200) in the for- 
mer . The interaction among the charge in the differ- 
ent channels produces a significant renormalization of the 
strength of the Coulomb interaction. This can be under- 
stood in the bosonization approach developed above, if 
we consider that each metallic nanotube in the rope con- 
tributes with four units to the number N of channels that 
can be bosonized jl9| . A large number N implies that the 
contribution of the factor C{x,t) to the supercurrent is 
greatly diminished, which has in practice the same effect 
as reducing the strength of the interaction. The picture 
is more involved considering the whole number of nan- 
otubes in the rope, as the interaction with the charge 
in the semiconducting tubules cannot be completely ne- 
glected. The overall physical effect can be taken into 
account by assuming a scale-dependent renormalization 
of the coupling constant 2e 2 /(f 2 vf) in the rope from its 
bare value at the scale of a few angstroms, as discussed 
in Ref. Q. 

Given that the total length of the sample RO3 is 
ps 1.7/zm, we have estimated the supercurrent by the de- 
cay of Tl through a distance L = 50/ k c pa 1.5/im. The 
evaluation of the prefactor in front of Eq. (^) is now 
more delicate, compared to that for the sample ST±. On 
the one hand, we have to bear in mind that the value of 
k c decreases according to the increase in the diameter of 
the sample. On the other hand, there are more metallic 
nanotubes in the sample RO3, in a number that may be 
estimated as 1/3 of the total number, which gives pa 60 
metallic nanotubes. Balancing both points, it is appro- 
priate to take now a prefactor in Eq. (|^) that is four 
times the value for the single- walled nanotube. 

We show in Fig. ^ the plots of Tl(T), including the re- 
sults obtained by adding the short-range attractive inter- 
action. The critical current in the sample RO3 at T pa OK 
is X pa 2.5/iA. We observe that the correct order of mag- 
nitude can be obtained from our results by considering 
a renormalization of the coupling 2e 2 /(tt 2 vf) down to a 
value pa 0.2, together with the effect of a weak short-range 
attractive interaction with coupling g/iirvp) ~ 0.15. As 
a final check of the consistency of our approach, we ob- 
serve that the curves in Fig. ^| reproduce the depen- 
dence on temperature measured experimentally in the 
sample -RO3, with the characteristic inflection point and 
the very slow decay around the critical temperature at 
T~10- 2 £U|. 
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of nanotubes M). In a sample like RO3, it can be al- 
ready observed that the supercurrent measured experi- 
mentally does not vanish near T c , which is at odds with 
the conventional picture of the proximity effect but in 
accordance with the results of our model. This enhance- 
ment of the superconducting correlations should deserve 
further study, in order to understand the experimental 
conditions under which the effect of the short-range at- 
traction may dominate over the Coulomb repulsion. 

Fruitful discussions with S. Bellucci and F. Guinea 
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supported by CICyT (Spain) and CAM (Madrid, Spain) 
through grants PB96/0875 and 07N/0045/98. 



FIG. 3. Plots of the critical current (in units of evFk c ) 
versus T/E c for different strengths of the interaction. 
From top to bottom, the solid curves correspond to 
2e 2 /(tt 2 vf) — 0.05,0.2, and the dotted (colored) lines to the 
latter interaction corrected by the additional short-range at- 
traction with coupling ig/^irvp) = 1.0,0.75,0.5 . 

We point out that the value 2e 2 / (tt 2 vp) » 0.2 is actu- 
ally quite close to what is predicted by integrating out the 
high-energy electron modes across two orders of magni- 
tude, from the angstrom scale to the diameter of the rope 
(w 20nm). This follows from the renormalization group 
approach worked out in Ref. Jf]J for the limit of a very 
large number N of subbands. Moreover, the required 
value of g matches what is expected from the coupling to 
the elastic modes of the nanotube. The short-range effec- 
tive attraction can be estimated from the modulation of 
the hopping t' = dt/da w 4.2eVA , the speed of sound 
v s s» 2.1 x 10 4 ms _1 , and the mass M of the atoms. This 
gives g/vp ~ t' 2 a 3 / (Mv 2 Vp) ~ 0.2, which is of the same 
order of magnitude needed in our fit. 

To summarize, we have seen that a model of interact- 
ing electrons in one dimension is able to explain the order 
of magnitude as well as the temperature dependence of 
the critical currents in both the ST\ and the RO3 samples 
of Ref. [p|- Our description is free of the shortcomings 
arising from the conventional picture of the proximity ef- 
fect, which relates the value of the critical supercurrent 
to the gap A and the normal resistance R through the 
expression 1 — ivA/(eR). Our approach focuses on the 
strong correlations in the one-dimensional electron sys- 
tem, explaining in this way why the experimental data of 
the supercurrent do not follow in general the temperature 
dependence of the gap in the superconducting contacts. 

Our discussion also stresses the relevance of the cou- 
pling to the elastic modes of the nanotube, which reveals 
itself through the presence of a short-range attractive 
electron interaction. This is also supported by recent 
experiments on the intrinsic superconductivity of ropes 



[1] J. W. G. Wildoer et al, Nature 391, 59 (1998). T. W. 

Odom et al, Nature 391, 62 (1998). 
[2] J. W. Mintmire, B. I. Dunlap and C. T. White, Phys. 

Rev. Lett. 68, 631 (1992). N. Hamada, S. Sawada and A. 

Oshiyama, Phys. Rev. Lett. 68, 1579 (1992). 
[3] M. Bockrath et al, Nature 397, 598 (1999). Z. Yao et 

al., Nature 402, 273 (1999). 
[4] J. Solyom, Adv. Phys. 28, 201 (1979). F. D. M. Haldane, 

J. Phys. C 14, 2585 (1981). 
[5] L. Balents and M. P. A. Fisher, Phys. Rev. B 55, R11973 

(1997). 

[6] Yu. A. Krotov, D.-H. Lee and S. G. Louie, Phys. Rev. 

Lett. 78, 4245 (1997). 
[7] R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 

(1997); Eur. Phys. J. B 3, 281 (1998). 
[8] H. Yoshioka and A. A. Odintsov, Phys. Rev. Lett. 82, 

374 (1999). 

[9] A. Yu. Kasumov et al, Science 284, 1508 (1999). 
[10] A. F. Morpurgo et al, Science 286, 263 (1999). 
[11] S. Bellucci and J. Gonzalez, Eur. Phys. J. B 18, 3 (2000). 
[12] J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Phys. 

Rev. B 59, 2474 (1999). 
[13] R. Egger and H. Grabert, Phys. Rev. Lett. 79, 3463 
(1997). 

[14] C. Kane, L. Balents and M. P. A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 

[15] D. L. Maslov, M. Stone, P. M. Goldbart and D. Loss, 
Phys. Rev. B 53, 1548 (1996). 

[16] R. Fazio, F. W. J. Hekking and A. A. Odintsov, Phys. 
Rev. B 53, 6653 (1996). 

[17] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 
1220 (1992); Phys. Rev. B 46, 7268 (1992); Phys. Rev. 
B 46, 15233 (1992). 

[18] The small relevance in general of intertube hopping in 
ropes has been discussed by A. A. Maarouf, C. L. Kane 
and E. J. Mele, Phys. Rev. B 61, 11156 (2000). 

[19] This approach has been used in the description of multi- 
walled nanotubes by R. Egger, Phys. Rev. Lett. 83, 5547 
(1999). 

[20; 



M. Kociak et al, report :ond-mat/0010220 



